All Maximally Entangled Four Qubits States 
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We find an operational interpretation for the 4-tangle as a type of residual entanglement, some- 
what similar to the interpretation of the 3-tangle. Using this remarkable interpretation, we are 
able to find the class of maximally entangled four-qubits states which is characterized by four real 
parameters. The states in the class are maximally entangled in the sense that their average bipartite 
entanglement with respect to all possible bi-partite cuts is maximal. We show that while all the 
states in the class maximize the average tangle, there are only few states in the class that maximize 
the average Tsillas or Renyi a-entropy of entanglement. Quite remarkably, we find that up to local 
unitaries, there exists two unique states, one maximizing the average a-Tsallis entropy of entangle- 
ment for all a > 2, while the other maximizing it for all < a < 2 (including the von-Neumann 
case of a = 1). Furthermore, among the maximally entangled four qubits states, there are only 3 
maximally entangled states that have the property that for 2, out of the 3 bipartite cuts consisting 
of 2-qubits verses 2-qubits, the entanglement is 2 ebits and for the remaining bipartite cut the en- 
tanglement between the two groups of two qubits is lebit. The unique 3 maximally entangled states 
are the 3 cluster states that are related by a swap operator. We also show that the cluster states are 
the only states (up to local unitaries) that maximize the average a-Renyi entropy of entanglement 
for all a > 2. 

PACS numbers: 03.67.Mn, 03.67.Hk, 03.65.Ud 



I. INTRODUCTION 

Entanglement lies at the heart of quantum physics. 
It was clear immediately after the discovery of quan- 
tum mechanics that entanglement is not "one but rather 
the characteristic trait of quantum mechanics, the one 
that enforces its entire departure from classical lines of 
thought" [![. Nevertheless, it was not until recently, that 
entanglement, besides of being interesting from a funda- 
mental point of view, was also recognized as a valuable 
resource for two-party communication tasks such as tele- 
portation [2[ and superdense coding [3]. With the emer- 
gence of quantum information science in recent years, 
much effort has been given to the study of bipartite en- 
tanglement [H; in particular, to its characterization, ma- 
nipulation and quantification It was realized that 
maximally entangled states are the most desirable re- 
sources for many quantum information processing (QIP) 
tasks. While two-party entanglement was very well stud- 
ied, entanglement in multi-party systems is far less under- 
stood, and even the identification of maximally entangled 
states in multi-party systems is a highly non-trivial task. 

The understanding of highly entangled multi-qubit 
states is crucial for the implementation of many QIP 
tasks in quantum networks. Highly entangled multi- 
qubit states, such as the cluster states or graph states, are 
the key resource of one-way or measurement based quan- 
tum computer Q, and as such raised enormous interest 
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in the QIP community. Even experimental realizations 
of one-way quantum computing with four-qubit cluster 
states has been demonstrated successfully Q. Highly 
entangled multi-qubit states are also the key ingredients 
of various quantum error correction codes and quantum 
communication protocols (§l-tlo|. However, unlike bipar- 
tite entanglement, very little is known about the charac- 
terization of entanglement in multi-qubits systems. 

The complexity in the characterization of entangle- 
ment in multi-partite systems can already be seen in 
the fact that for 3 qubits there are essentially two types 
of genuine 3-partite entanglement and even the notion 
of maximally entangled states is not unique [ill ]. One 
can think of both the GHZ-state, |GHZ) = (|000) + 
|111))/V2, and the W-state \W) = (|001) + |010) + 
|100))/v3 as two types of maximally entangled states, 
that are not related by stochastic local operations and 
classical communication (SLOCC). Nevertheless, in the 
case of 3-qubits, one can single out the GHZ-state as 
the unique maximally entangled state for the following 
two reasons. First, the GHZ class (i.e. the set of states 
that can be obtained from the GHZ state by SLOCC) 
is dense in the space of 3-qubits and therefore the W- 
class is of measure zero. This means that by LOCC it is 
possible to convert the GHZ state to a state that is arbi- 
trarily close to the W-state, but the W-state can not be 
converted (not even by SLOCC) to a state that is close 
to the GHZ-state. Second, the GHZ-state is the only 
3-qubit state with the property that the bipartite entan- 
glement between any one qubit and the other two qubits 
is maximal; that is, the reduced density matrix obtained 
after the tracing out of any two qubits is proportional to 



2 



the identity. 

Similarly, for n-qubits one can define maximally entan- 
gled states as states with the property that the reduced 
density matrix obtained after the tracing out of any k 
qubits, with n/2 < k < n — 1, is proportional to the iden- 
tity. For example, the codeword states of the 5-qubits er- 
ror correcting codes are maximally entangled 12]. How- 
ever, as we show below, for four qubits such states do 
not exist. It is also known [l3| that maximally entangled 
states exist for n = 6 and do not exist for n > 8. To the 
authors knowledge, the case of n — 7 is unknown. 

In this paper we find an operational interpretation of 
the 4-tangle which enable us to characterize all maxi- 
mally entangled four-qubits states. We define a state 
to be maximally entangled if its average bipartite en- 
tanglement with respect to all possible bi-partite cuts is 
maximal (e.g. see [13, EBl and references therein). More 
precisely, we divide the four qubits into two groups, each 
consisting of two qubits, and calculate the pure bipar- 
tite entanglement between the two groups of qubits. We 
then find the class of all states that maximize the aver- 
age entanglement of the 3 = ( 2 )/2 such (inequivalent) 
bi-partite cuts. We find that when we take the mea- 
sure of bipartite entanglement to be the tangle, there is 
a 4-real parameter class of states M. that maximize the 
average tangle. However, when we take the measure to 
be the entropy of entanglement, or the Tsallis and Renyi 
a— entropy of entanglement, we get that up to local uni- 
tary there are only two states that maximize the average 
a— entropy of entanglement. Quite remarkably, we find 
that up to local unitaries the state 

\L) = [n + wui + uj 2 u 2 ] , (1) 
where uj = e l27r / 3 , 

u = \4> + )\4>+) , Ul = |0-)|0-) 
u 2 = , u 3 = \ip-)\ip-} 

and = (|00)±|ll))/V2and 1^) = (|01>±|10))/V2, 
is the only state that maximize the average Tsallis 
a— entropy of entanglement for all a > 2 (see Fig. 1). In- 
teresting properties of the state \L) have been discussed 

in |il Eg. 

On the other hand, we show that the state 

|M) " 71 Uo + 7! [ui + U2 + M3] (2) 

is the only state that maximize the Tsallis a— entropy of 
entanglement for all < a < 2. Ten years ago the state 
\M) was conjectured to maximize the entropy of entan- 
glement [17|. More recently, it was proved that locally 
it is indeed maximally entangled [l8j |. In Fig. 1 we draw 
a graph of the average Tsallis a-entropy of entanglement 
as a function of a for the states \M) and \L). 

In addition, among the maximally entangled four 
qubits states, we identify 3 ultimate maximally entan- 
gled states that have the property that for 2, out of the 
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FIG. 1: A graph of the average Tsallis a-entropy of entan- 
glement as a function of a. The blue line corresponds to the 
state \M), the green line to the state \L), and the dashed red 
line to the cluster states. Like the cluster states, the graph 
for any maximally entangled state in M is between the blue 
and green lines. 



3 bipartite cuts, the entanglement is 2 ebits and for the 
last bipartite cut the entanglement between the groups 
of two qubits is lebit. The unique 3 maximally entangled 
states are the 3 cluster states that are related by a swap 
operator (but not by SLOCC): 



Id) = - (10000) + 
\C 2 ) = \ (10000) + 
\C 3 ) = \ (10000) + 



|1100> 4 


- 10011) 


-11111)) 


(3) 


|0110) 4 


- 11001) 


-11111)) 


(4) 


|1010) 4 


- 10101) 


-11111)) 


(5) 



We show that these cluster states are the only states that 
maximize the Renyi a— entropy of entanglement for all 
a > 2. 

This paper is organized as follows, in section (II) we 
discuss the generic class of four qubits states, consist- 
ing of an uncountable number of SLOCC-inequivalent 
classes. In section (III) we find an operational interpre- 
tation of the 4-tangle and discover a four real parameter 
class of all four qubits states that maximize the aver- 
age tangle. We then use this result in section (IV) to 
find maximally entangled states with respect to other 
measures of entanglement, such as the Tsallis and Renyi 
a— entropy of entanglement. In section (V) we discuss 
more maximally entangled four qubits states. We end 
in section (VI) with a summary, conclusions and a dis- 
cussion on the extension of the results presented here in 
higher dimensions. 
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II. UNCOUNTABLE NUMBER OF FOUR 
QUBITS SLOCC-INEQUIVALENT CLASSES 

In [2l| it was argued that 4-qubits pure states can be 
classified into nine groups of states. One of these nine 
groups is called the generic class as with the action of 
SLOCC it is dense in the space of 4-qubits H4 = C 2 ® 
C 2 <g> C 2 ® C 2 . The generic class is given by 



A 



Z U + 011*1 + z 2 u 2 + z 3 u 3 



z 0l z 1 ,z 2 ,z 3 G Cj . 



In [2(j, [21| it has been shown that all the states that 
are connected to the class A by SLOCC form a dense 
set of states. That is, the class of states GA, where 
G = SL(2, C) ® SL(2, C) <g> SL(2, C) ® SL(2, C), is dense 
in ^4. In the following we discuss several properties of 
the generic class that will be very useful for our theorems 
in the next sections. 

For k = 0, 1, 2, 3, we denote by \k)) = \ij), with i, j — 
0,1, a state of two qubits, such that ij is the binary 
representation of k. Hence, any \ip) G %4 can be written 
as 



(6) 



k=0 fc'=0 



where {|fc))} is the computational basis of qubits 1 and 2, 
and is the computational basis of qubits 3 and 4. 

With these notations we define the following four quan- 
tities: 



Definition 1. Let 



G Ua. Then, 



Det [Tip] 



Tr 



/ \ ? 



if 771 = 

if 777 = 1.2,3 



where is the 4x4 matrix whose components T^w are 
defined in Eq. ©, and 



J = 



1 
0-10 
0-100 
10 



The four polynomials defined above take a simple form 



on A. If 



then 



zo ziz 2 z 3 



2m 




Z 



2m 



+ z\ m + Z- 



2)n 



if 777 = 

if m = 1,2,3 



In 20, 21] it has be shown that these four pohnomials are 
invariant under the action of the group G = SL(2, C) ® 
SL(2, C)®SL(2, C)®SL(2, C). That is, if g € G and \tp) G 
then £ m {g\fl) = £m{W)) for all 777 = 0, 1, 2, 3. Other 
polynomials that corresponds to true 'tangles' have been 
considered for example in [lj| . As discussed in [U |22j , 
one of the consequences of this property is that the four 



functions f m = \E V 



(777 = 1,2,3) and / 



entanglement monotones 34|. However, here we show 
that this property implies that almost all the states in A 
are not related by SLOCC, which means that A contains 
an uncountable number of SLOCC inequivalent classes 
of states. 

Proposition 1. (c.f. fM, HI Let \ip) and \ip') be 
two normalized states in A. Then, the transformation 
\ip) — !• \ip'} can be achieved by SLOCC only if f m — 
f m (|^)) for allm = 0,1, 2, 3. 

Proof. Let ip G A be a normalized state and let g G G. 
Then, 



fn 



\g^\\ 



fm(g^) fm(lp) 



\91p\\ 



Ifl'V'll 



where in the last inequality we used the fact that for 
ip G -4 and g e G, \\gip\\ >M>\\ = 1 (see Appendix E] for 
the Kempf-Ness theorem [2J|). □ 

In the definition above E m is defined only for m < 
3. The absolute value of the polynomials with higher 
values of m are also entanglement monotones, but they 
are in the algebra generated by these four polynomials 
and therefore do not contain any additional information 
about the entanglement of the states. 

The proposition above implies that the class A con- 
tains an uncountable number of states that are not con- 
nected by SLOCC transformation. More precisely, if 
\ip) = J2j=o z j u j an d W) = Y^j=o z 'j u j^ and there is 
no permutation a such that Zj = ±2^.^ for j = 0, 1,2, 3 
with an even number of — signs, then the transforma- 
tion j^) — > can not be achieved by SLOCC. More- 
over, from the proposition below it follows that on A if 
two states are connected by SLOCC operation, then the 
transformation must be a local unitary. 

Proposition 2. fM, \M, \M] Let \tp') G A. Then, 
the transformation — > can be achieved by a local 
unitary U G SU(2) <g> SU(2) <g> SU{2) <g> SU{2) if and only 

if£ m {m=£m{W))- 

For the purpose of this work, we generalize the propo- 
sition above to include all local unitaries; that is, not only 
those in SU(2) (g) SU(2) ® SU(2) ® SU(2). 



Proposition 3. Set 



fi — l^i — £0 



h = \el~e 2 \\ h = \sl-s 3 \ 2 



Let \ip), \ip') G A. Then, the transformation \tp) \tp') 
can be achieved by a local unitary U G U(2) ® U(2) i8> 
U[2) ® U(2) if and only if f m (|t/>)) = f m for all 

integers < 777 < 6. 

Proof. Note that the first 4 conditions (i.e. m = 0, 1, 2, 3) 
imply that 

£ (VO = oW) . £iW = b£i(il/) 
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with |o| = \b\ = |c| = \d\ = 1. / 4 (» = fi{ip') implies that 
a = b 2 . Similarly, the condition on / 5 implies c = 6 2 , and 
the condition on /g implies d — b 3 . Now, write 6 = r 2 . 
We therefore have 

£oW>) = SoW) , £ 1 (^) = £ 1 (r^') 

Thus, from proposition^ \tp) and r\ij/) are related by 
a local unitary in SU(2) ® SU(2) ® SU(2) ® SU(2). The 
argument clearly can run backwards. □ 

As we show now, among the 4 entanglement monotones 
f m (m = 0, 1, 2, 3), the 4-qubits entanglement monotone 
/1 is the only one that is invariant under any permutation 
of the 4-qubits. In fact, we find that this monotone is the 
4-tangle. 

A. The monotone /1 = |fi| and the 4-tangle 

Given a bipartite state \ip AB ) € C" <g) C m , the tangle 
(or the square of the I-concurrence) is defined by 

t ab = r (|<0) = S L = 2 (1 - Trp 2 ) , (7) 

where p r — TrB\tp AB ) (tp AB \ is t ne reduced density ma- 
trix, and Sl is the linear entropy. 

For two qubits the tangle can be expressed as the 
square of the concurrence; that is, for \ip) 6 C 2 <£> C 2 
the tangle is 

tab = |(¥#)| 2 where \ip) = cr y (g> a y \ip*) , 

and a y is the second Pauli matrix. Note that the basis 

is chosen such that o~ y — ( ^ q ~\ , and if in this basis 

\f) = T,i,j a ij\ij) then \<P*) = Ei,i a *iNi)- 

For mixed two qubits state p AB , the tangle is defined 
in terms of the convex roof extension: 

tab = t(p ab ) = min^pjr(l^i)) , 

i 

where the minimum is taken over all the decompositions 
of the form p AB = Y,iPi\i>i){i } i\- 

In [26(1 it was shown that one can extend the definition 
of the 2-qubits tangle to 3-qubits. Given a 3-qubits pure 
state \ip) £ C 2 ® C 2 (8 C 2 the 3-tangle is defined by 

TABC = T (\lp)) = T A (BC) - TAB - TAG , 

where tab = t(p ab ) (with p AB = Trcl^iV 1 !): an d 
ta(bc) is the tangle between the qubit system A and 
two qubits system BC. In (2(| it was shown that the 3- 
tangle is non-negative and its square root has been proved 
to be an entanglement monotone in (22|. It was also 
shown [26] that it is symmetric under permutations of 
the three qubits A, B, and C. From its definition, the 
3-tangle can be interpreted as the residual entanglement 



between A and BC, that can not be accounted for by the 
entanglements of A and B, and A and C, separately. 

The Wong-Christensen 4-tangle [27j is defined simi- 
larly. Let € Ha = C 2 <g> C 2 (8) C 2 <g> C 2 , the 4-tangle is 
defined by (27| 

TABCD = \ (i>\<Ty <8) <T y <Z> <Jy <8) (T y \lp*)\ 2 . 

In [13] the 4-tangle was shown to be an entanglement 
monotone and invariant under permutations. In the next 
section we will see that like the 3-tangle, the above 4- 
tanglc can also be interpreted as a type of residual en- 
tanglement. Moreover, as we show now, the square of 
the monotone /1 is the Wong-Christensen 4-tangle. 

Proposition 4. Let 6 H4. Then, 

TABcoim - mm 2 ■ 

The above proposition follows directly from the fact 
that there is a single SL(2 } C)® 4 invariant polynomial 
with homogeneous degree 2 [28j . Since both the 4-tangle 
and \£± | 2 have these properties they must be equal. In 
the proof below we show this equivalence by a direct cal- 
culation. 

Proof. Denote \ip) = ^/po"|O)|0°) + ^I)^ 1 ) , where 
€ C 2 ® C 2 ® C 2 with j = 0, 1 are three qubits or- 
thonormal states. Also, denote 

w)= E a * (J) i i > 

ie{o,i} 3 

for j = 0, 1. With this notations it is straightforward to 
show that both £i(\ip)) and equals 

9 fprrrfJ K^ ^n^h,^ +n (0 hM*> j-J°K^ 

^VPOPl \^ U 000 U 111 ' a 110"001 "1 a 101 a 010 u 011 u 100 
"lll u 0Q0 u 001 u 110 u 010 u 101 u 100 u 01lJ- \°) 

This completes the proof. □ 

III. OPTIMIZING THE AVERAGE TANGLE 

As a measure for pure bipartite entanglement we first 
take the tangle or the square of the I-concurrence (see 
Eq-©)- Now, in four qubits there are 4 bipartite cuts 
consisting of one-qubit verses the rest 3 quibts and 3 bi- 
partite cuts consisting of 2-qubits verses the rest 2 qubits. 
Denoting the four qubits by A, B, C, and D, we define 

T\ = 7 ( T A(BCD) + T B(ACD) + T C(ABD) + T D(ABC)) (9) 
Tl = - (T(AB)(CD) + T{AC)(BD) + T(AD)(BC)) ) ( 10 ) 

where ta(bcd)i f° r example, is the tangle between qubit 
A and qubits B,C,D. Similarly, t/ab)(gd)> f° r example, is 
the tangle between qubits A,B and qubits C,D. Note that 
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the maximum possible value for n is 1 and the maximum 
possible value for T2 is 3/2 (since a maximally entangled 
4x4 bipartite state has tangle 3/2). However, it is argued 
now that no 4 qubit pure state can achieve this value for 

In fijl, it has been shown that 



n<T2< -ti 



(11) 



Hence, since t\ is bounded by 1, it follows that r 2 < 
4/3 < 3/2. That is, there are no 4-qubit states for which 
all the 3 reduced density matrices, obtained by tracing 
out two qubits, are proportional to the identity. More- 
over, from the inequality above, it follows that for states 
with r 2 = 4/3, T\ must be equal to 1. In the following 
theorem, we characterize all states with t\ = \. 

Theorem 5. fM,\M] Let \ ip) £ H4 = C 2 <g> C 2 <g) C 2 <g) C 2 
be a normalized four qubit state. Then, 

r i (IV')) = 1 */ an d on l y tf G -4. , 
up to local unitary transformation. 

A weaker version of the theorem above has been first 



pointed out in [2l|. In [2l| the authors argued that 



among all the states in GA, only states in A have t\ = 1. 
A year later in [22[ theorem [5] was fully proved. Never- 
theless, for the purpose of completeness, we provide here 
a proof of TheoremJS] for all states in ^4, independently 
of the work in [2lL|22|. 

Proof. Using the Kempf-Ness theorem [24| (see also Ap- 
pendix [A} applied to G (defined above) and the fact that 
GA is the set of stable vectors 35[, one can show that 
ip £ KA, with K = SU{2) <g> SU{2) <g> SU{2) ® SU(2), if 
and only if [U 



(tP\X\tP)=0, 



(12) 



for all X in Lie(G). Note that Lie(G) acting on W4 is 
the direct sum of the Lie algebras of SL(2,C) acting on 
one tensor factor. Now, let ip be a normalized state with 
T~i(ip) — 1. Therefore, we can write ip = \0)\<po) + \l)\<pi), 
with (cpo \ipo) =1/2 = (ip\\<pi) and (volvi) = 0. Low, let 



U 



a b 

-b* a* 



be a unitary matrix with \a\ 2 + \b\ 2 = 1. Then, if U\ is 
U acting only in the first tensor factor we have 



Thus, the maximum value is attained at U = I. This 
implies that the condition in Eq.([T2|) is true for the part 
of the Lie algebra coming from the elements that act only 
on the first factor. The argument for the other factors is 
the same. □ 



From the Eq. (ITT1) it follows that among all the states 
with T\ = 1, we have 

4 

1< T 2 < - . 

- - 3 

It is interesting to note that the 4-qubit GHZ state gives 
the minimum possible value for r 2 . That is, it is the 
least entangled state among all the states with n = 1. 
On the other hand, for all the 3 cluster states defined 
above, r 2 = 4/3. In fact, as we will see later, the cluster 
states are the only states that achieve the maximal value 
for r 2 in such a way that 2 of the terms (i.e. tangles) 
appearing in the definition of r 2 (see Eq. (fTU)) ') are equal 
to 3/2 and one of the terms equals to 1. 

From Theorem[5]and the inequality (fTTj). it follows that 
only states in A can maximize r 2 . From the following 
theorem it also follows that states that maximize t 2 must 
have zero 4-tangle. 

Theorem 6. Let tp £ be a 4-gubits pure state and 
denote by tabcd{^) its ^-tangle (defined above). Then, 



T2«0 



4ti(V>) - t~abcd(iP) 



Remark. The equation above can be written as tabc d — 
4ti — 3r 2 , where At\ can be interpreted as the total 
amount of entanglement in the system, whereas 3t 2 
can be interpreted as the total amount of entanglement 
shared among groups consisting of two qubits each. In 
this sense, the 4-tangle can be interpreted as the resid- 
ual entanglement that can not be shared among the two 
qubits groups. Note that from the equation above it is ob- 
vious that the 4-tangle is invariant under permutations. 

Proof. Following the same notations as in Proposition [4j 
we denote \i\>) = y/po\0)\(j) ) + y/pill)]^ 1 ) , where the 3- 
qubits states (qubits BCD), \<fp) £ C 2 <g> C 2 <g> C 2 with 
j = 0,1, are orthonormal. We also denote by D x with 
X £ {B, C, D} the discriminant of \ip) [l9l |: 



where 



af = Tr^ x |0 fc )(0 fc ' I with k, k' £ {0, 1} 



X> x = Tr (er^CTx 1 



7 x a x) 



and the trace is taken over all the remaining two qubits 
that are not the X qubit. The sum of the discriminants 
is denoted by D = D B + T>c + ^>d- 

With these notations the one-qubit reduced density 
matrices can be written as follows: 



/ = Tr #A |^| =Po|0)(0| 
p X = Tr^)(V| =P0<Tx+Pi° 



(13) 



Similarly, the two-qubit reduced density matrices are 
given by 



=Tr* AX \1>){4>\= J2 VPkW\k){k'\®o- k x k ' . 
fc,fc'e{o,i} 

(14) 
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Substituting these reduced densities matrices in the ex- 
pressions for the linear entropy gives 



Since Sl(p a ) — ^PoPi: summing over X gives 

4r 1 -3r 2 = 4poFi(l-S) • 
Now, if we denote 



E 

ie{o,i} 3 



aV\i) 



for k = 0, 1 , then a straightforward calculation (sec 
Eq.(3.38) in 0) gives 



D = 1 



_(o) (i) , Jo) (i) 

u 000"lll ' a 110 u 001 



(0) (i) 
a 101 a 010 



(0) (1) 
a 011 a 100 



(0) (1) 

a m a ooo 



(0) (1) 
a 001 a 110 



„(0) (1) 
a 010 a 101 



(0) (1) 
a 100 a 011 



(15) 



A comparison of this expression with the one given in 
Eq.© implies that 4ti — 3r 2 = tabcd- □ 

From Theorem [5] and Theorem [6l we have the following 
corollary. 

Corollary 7. A normalized state \ip) € 7H4 is maximally 
entangled (i.e. T2(\ip)) = 4/3J if and only if up to local 
unitary \ip) € M., where M. is the set of states in A with 
zero 4-tangle. 



A state ijj = 
parameters Zj (j 



ZjUj in A depends on four complex 
0,1,2,3). The condition that the 

4-tangle tabcd (VO = IXX-=o 2 J| 2 = implies that the 
states in the maximally entangled class M. are character- 
ized by 4 real parameters since we also have the normal- 
ization condition and we ignore the global phase. If we 

jo 



write Zj = y/p]e 1 j in its polar form (with non-negative 
Pj and 9j € [0, 2ir}) we can characterize the class A4 as 
follows: 



M 



/pjt 



3 

Eft 

3=0 



1 1 E p j ( 

3=0 



(16) 



Minimization of ti 



From theorem [SJ it also follows that the states in A 
with the minimum possible value T2 = 1, can be charac- 
terized as follows. Denote by T m in the class of all such 
states. Then, 



7I„in = {^e^|r 2 (V') = l} 



3 3 

E^' u j| E x j = 1 ' x i e 

3=0 3=0 



(17) 



Note that the four qubits GHZ state belongs to Tmin- In 
this sense, the GHZ state is a state in A with the least 
amount of entanglement. 



IV. DIFFERENT MEASURES OF 
ENTANGLEMENT 



Up to now we took the measure in Eqs. (|9ll0p to be 
the tangle, which is given in terms of the linear entropy. 
The measures of entanglement that we consider in this 
section are Renyi entropy of entanglement and Tsallis 
entropy of entanglement. We denote these measures by 
and , respectively. These bipartite measures of 
entanglement are defined as follows. Given a bipartite 
state \ip AB ) G C" <g> C m , the entanglements E^ and 
E^\ measured by the Renyi and Tsallis entropies, are 



given by 



R 



(a) 



{\^ AB )) = - logTrp? 

1 — a 



E^ (|VO) = ^- (Trp« - 1) 

where p r = TiB\ip AB )(iJ AB \ is the reduced density ma- 
trix and the log is base 2. Note that both Renyi and 
Tsallis entropies approach the von-Neumann entropy in 
the limit a — > 1. The Tsallis entropy is concave (see 

for example [2!|) and therefore E^ is an ensemble en- 
tanglement monotone (i.e. non-increasing on average 
under LOCC). The Renyi entropy is also concave for 
< a < 1, but only Shur concave for a > 1 [3(| ■ Hence, 

for a > 1, E^ is only a deterministic monotone (i.e. 
non- increasing under deterministic LOCC). Nevertheless, 
unlike the Tsallis a-entropy of entanglement with a^l, 
the Renyi a-entropy of entanglement is normalized nicely 
so that it is equal log d for maximally entangled states in 



Similar to the definition of the average tangles (|9|10|) . 
we define the average a— entropy of entanglement of four 
qubits states as follows: 



E 



(a) 



E. 



(a) 



i( 
U 



where E^ BCD y for example, is the Tsallis or Renyi a- 
cntropy of entanglement between qubit A and qubits 
B,C,D, where it will be clear from the context if we mean 
Renyi or Tsallis. Note that due to Eq.(fTTj) the maximum 

value of E% cannot be 2 ebits (i.e. the same value as 
the the value for two bell states). 

From Theorem [5] we know that E[ a '(\ijj)) = 1 iff G 
A. However, for the general measures of entanglement we 
do not have an equation analog to Eq. (fTTj) and therefore, 
can not argue that if \ip) maximize E^ then it must 

maximize E[ a '. Nevertheless, we will see that for the 
average a— entropy of entanglement with a > 2 this is 
indeed the case 1361] - 



,(a) 

A(BCD) 



(a) 

(AB)(CD) 



E 



(a) 

B(ACD) 



E, 



(a) 

C(ABD) 



D(ABC) J 

(18) 



E 



(a) 

(AC)(BD) 



E 



(a) 

(AD)(BC) 



(19) 
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In order to optimize E 2 a ^ , we first prove the following 
theorem. 

Theorem 8. Let p be a 4 x 4 normalized density matrix 
with eigenvalues Ao > Ai > A2 > A3. Let Sl{p) = t 

and denote x(t) = ^Jl — |r and y(r) = yjl — |r. The 

maximum (minimum) possible value of either Tsallis or 
Renyi entropies with < a < 2 (a > 2) is obtained if 
and only if is given by 



1 + 3x(t) l-x(r) l-x(r) I — x(t) 



(20) 



The minimum (maximum) possible value of Tsallis or 
Renyi entropies with < a < 2 (a > 2) is obtained if 
and only if the set {Xi} is given by 



1 + x(t) 1 + x(t) 1 + x{t) 1-3o;(t)1 4 3 
4 ' 4 ' 4 ' 4 I ^3- T -2> 



1 + V{r) 1 + yjr) 1 - 2y{r) 
3 ' 3 3 



(I < T < -) 

3 

(0 < T < 1) 

(21) 



One direction of theorem [5] has been proven in [3l[. To 
complete the proof of theorem [51 we show in Appendix 
|B1 that the Renyi and Tsillsa entropies obtain there ex- 
ternum values only for the sets of eigenvalues that ap- 
pear in the theorem. 



qubit normalized state £ H4, with T( 12 )(34) = ti, 



F (13)(24) 



= t2, and T| 



(14) (34) 



f^lihMM) < E^\m < f£l(h,t 2 ,t 3 ) . (22) 

From theorem [5] and the definition of Tsallis entropy it 
follows that for a > 2 



^max(^) — 
1 



1 



■ X 



a — 1 

3(l+a:) n _ (l-3a:)° 



1 - 



(i-M! 



for § < t < § 



for 1 < t < § 



(23) 



where x — x(t) = ^Jl — |t and y ~ y(t) = ^ 
Furthermore, from theorem [5] it follows that for < a < 
2, E^ n (t) is given by the exact same expression as in 



Eq.(E2 



Note that Emax(t) and E^'^t) are continuous at the 
point t = 4/3 (i.e. x = 1/3 and y — 0), but not their 

derivatives. Nevertheless, since the derivatives of -Emax(i) 
and E^) n (t) in both regions 1 < t < 4/3 and 3/4 < t < 

3/2 are positive, it follows that Ema*(t) and E^ n (t) are 
both monotonically increasing with t. Now, from Eq. (fTTj) 
it follows that i 1+^2+^3 < 4 and therefore the RHS of the 
Eq. (1221) reach its maximum value when t\ + t 2 + t 3 = 4. 

Hence, for any -0 £ H4, E^ {ip) is bounded above by 
{a > 2) 



A. The state \L) 

In this section we show that the state \L) in Eq. (fT]) 
has the remarkable property that it maximizes the aver- 
age Tsallis entropy of entanglement E 2 a ^ for all a > 2. 
In addition, among all the states in M., the state \L) 
minimizes E 2 a ^ for all < a < 2. 

Theorem 9. (a) Let ^6^4. Then, 

E ( 2 a \ip) < E { 2 a \\L)) for alla>2 

with equality if and only if up to local unitaries if) — \L) . 
(b) LettpeM. Then, 

E ( 2 a \^) > E { 2 a \\L)) for all0<a<2 

with equality if and only if up to local unitaries if) = \L) . 

Proof. Let E^{t) and E^ n {t) be, respectively, 
the maximum and minimum values E^"\^j) 
can take among all normalized bipartite states 
-0 € C 4 <g> C 4 with tangle r(ip) = t. Fur- 



UM = max {/<&(*!, t 2 ,t 3 ) I ^i fc = 4 , t k < 3/2} 



fe=i 



Similarly, for any ip £ M., E 2 (tp) is bounded below by 
(0 < a < 2) 



3 

L<«> = min ta.tg) | £)f fc = 4, i fc < 3/2} 



fc=i 



ther, let finL(ti,t 2 ,t 3 ) = 1/3 Efe=l E )£L{t k ) and a < 2; ^ q) (|V>)) = £ (q) if and only if {Pj} = {Qj} 



In Appendix Owe show that = /iSL(4/3, 4/3, 4/3) 
and = 7^(4/3,4/3,4/3). We also show that the 
point ti = t 2 = t 3 = 4/3 is the only point of global max 
for fm&x and the only point of global min for fj^ n - 

Let £ V.4, and denote by {Pj}, {Qj}, and {Rj} 
the eigenvalues of the reduced density matrices of \tp) 
obtained after tracing out qubits C and D, B and D, 
and B and C, respectively. From the analysis above and 
the results in appendix [Cl E 2 a \\->p)) — if and only 
if T (i2)(34)(V0 = T(i2)(34)(V0 = T(i2)(34)(-0) = 4/3 and 
the distributions {Pj}, {Qj}, and {Rj} are all of the 
form given in Eq. (|2"Tj) with r = 4/3. That is, for a > 2 

E { 2 a) (\^)) = UW if and only if {P 3 } = {Qj} = {Rj} = 
{1/3,1/3,1/3,0}. Similarly, if £ M then, for < 



f£L(t u t 2i t 3 ) ^ l/3ELl^it(*fe)- Th us, for a 4 



{Rj} = {1/3,1/3,1/3,0}. 
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A priori it is not clear that such a state with {Pj} = 
{Qj} = {Rj} = {1/3,1/3,1/3,0} exists in % 4 . How- 
ever, now we show that up to local unitaries there exists 
exactly one state with this property and the state is \L). 

First note that if {Pj} = {Qj} = {Rj} = 
{1/3,1/3,1/3,0} then up to local unitaries \ip) € M. C 
A. Therefore, we can write \ip) = zqUq + ZiUi + Z2U2 + 
Z3W3. In this case, the eigenvalues of the reduced density 
matrices are given by (j = 0, 1, 2, 3): 



P 



Zk 



k=0 



R i 



k=0 



kZk 



(24) 

where Ajf. and Bjk are the matrix element of the two 
orthogonal 4x4 orthogonal matrices: 




Now, since {Pj} = {1/3, 1/3, 1/3, 0} we have 



n 1 
z Q = , z x = -j=e 



Z-2 



Z3 



were we have used the fact that the transformation \tp) — > 
\tp') = Ylj=o z <?(j) u i can t> e achieved by local unitaries for 



Ro 
Ri 
R2 



+ e 1 " 2 + e 
+ e 102 + e 



all permutations a. With these values for Zj we get 
1 

12 
1 

12 
1 



Q 3 = 



12 1 

6 12 1 



(25) 



Thus, {Qj} = {Rj} = {1/3,1/3,1/3,0} if and only if 
up to permutation e idl = ±1, e i02 = ±e"/ 3 , and e iH = 
_l_ e i27r/3 _ p rom proposition [5] it follows that up to local 
unitaries \ip) = \L). □ 



B. The State |M) 

In this section we show that the state \M) in Eq. @ 
has the remarkable property that it maximizes the aver- 
age Tsallis entropy of entanglement E% for all < a < 
2. In addition, among all the states in Ai, the state \M) 
minimizes E2 f°r all a > 2. 

Theorem 10. (a) Let ip € A. Then, 

E ( 2 a \iP) < E ( 2 a \\M)) for0<a<2 (26) 

with equality if and only if up to local unitaries ip = \M) . 
(b) Let ip G M. Then, 



E { 2 a \iP) > E ( 2 a \\M)) for a > 2 



(27) 



with equality if and only if up to local unitaries tp = \M) . 

Proof. Let ip € A, and denote Zj = y/pje l0j for j — 
0, 1, 2, 3. With this notations, £ x (ij}) = J2 j= o Pj ei2&1 > and 
for a fixed value of \£%\ = a, (a > 0), the formula in 
Theorem [5] can be written as 

t 1 (ip)+t 2 (ip)+t 3 (ip) = 4-a 2 . 

Moreover, note that if = a then pj < (l + a)/2 for 

all j = 0, 1, 2, 3. Therefore, we denote by E^x{a, t), the 
maximum value £lp (ip) can take among all normalized 



Z 4 with tangle t(<p) 
o)/2. 



t and 



bipartite states <p € C 4 yy ^ 
Schmidt coefficients pj < (1 

Now, a simple calculation shows that for a distribution 
of the formpo > Pi = P2 = P3, we get that pj < (l + a)/2 
if and only if t > 2(1 — a)(2 + a)/3. Therefore, from 
theorem [8] it follows that this is the optimal distribution 
for t > 2(l-o)(2 + o)/3. On the other hand, if t < 2(1- 
a)(2 + a)/3, the optimal distribution (up to permutation) 
is given by po = list > p\ > P2 = P3 . This is follows from 
the extension of the results in f3lf. and in particular [32| , 
it is a consequence of Eq.(22) in [3lj | . Hence, we conclude 
that for < a < 2 



E (a) (a t) = - 

a — 1 



3(1-t)° _ (l+3z)° 



for i > 



2(l-q)(2+q) 



1 _ I l+Q 'l" _ (l-°+2u)° _ 2(l-q-qj)° 
1 V 2 / 6° 6" 



otherwise 
(28) 



where x = x(t) = yl — |< and u = uj(t) = 

V4 - 2a - 2a? - 3t. 

Similar to the definition in theorem [9j we define 

f^{a,t x ,t 2 ,t 3 ) = 

i (Sil(a, tx) + EtL(a, t 2 ) + E^(a, t 3 )) . (29) 

We therefore have for < a < 2 

E ( 2 a) W < f^{aM,t 2 ,h). 

A straightforward calculation, similar to the one given 
in Appendix [Cj shows that the global maximum of the 

function /ma,x(a, ^1,^2, t 3 ) is unique and is obtained at the 
point a = (i.e. ip £ M) and t\ = t 2 = t 3 = 4/3. There- 
fore, from theorem [5] (and in particular from Eq. (|20[) 
with r = 4/3), it follows that this global maximum is 
obtained if and only if 

Pj, Qj, JJ,. G {iiii} Vj =0,1,2, 3, (30) 

where the sets {Pj}, {Qj}, and {Rj} are the eigenval- 
ues of the reduced density matrices of ip obtained after 
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tracing out qubits C and D, B and D, and B and C, 
respectively. A priori, it is not clear that a four qubits 
state with these properties exists. We show now that 
there exists only one state (up to local unitaries) with 
these properties, and it is the state \M). Note also that 
if such a ijj exists then ip £ M. C A. 

Up to a local unitary (see proposition [3]), w.l.o.g. we 
can assume that Pq = 1/2 and P\ = P2 = P's = 1/6. 
Further, due to the freedom of global phase, we have 

Next, the condition that three of the QjS and three of 
the RjS equal to 1/6 implies that \ip) equals to one of the 
8 states: 

T2 Ua± Ts Ul± Tz U2± Ts U3 

Using proposition |3] we get that all these eight states are 
equivalent under local unitaries. We are therefore left 
with one state 

™ 7T° + 7T 1 + 7T 2 + 7T 3 (31) 

Up to a local unitary, this state is the same as \M) in 
Eq. ©. 

The proof of part (b) of the theorem follows the same 
lines as above with a = and a > 2. □ 



The eigenvalues of the reduced density matrices are 
given by 0' = 0,1,2,3): 





3 


2 


3 






, Rj = 


BjkZk 




fc=0 




fc=0 



Hence, the equality H ({Pi}) — 2 leads to 
1 , 



2 6 



3=0, 1,2,3 



Now, since a quantum state is defined up to a global 
phase, w.l.o.g. we can take 9q = —6\ which reduces 
the number of free parameters to 3. Next, the condi- 
tion H({Qj}) = 2 together with proposition 2 implies 
that up to a local unitary 

1 1 e %1 e* 7 
\4>) = 2 U ° ~ 2 Ul + ~T U2 + T U3 ' 

where 7 is the only free real parameter left. Now, it is 
a simple calculation to check that at least two of the Rj 
(j = 0,1,2,3) equals to zero. Therefore, H ({Rj}) < 1. 
If we require that H ({Rj}) — 1, then e* 7 = ±i and up 
to a local unitary = IC2). □ 

Next we show that the three cluster states are the only 
states (up to local unitaries) that maximize the average 
Renyi entropy of degree a > 2. 



The cluster states 



1. Cluster states are the only states that maximize the 
average Renyi entropy with a > 2 



We are now ready to analyze the maximally en- 
tangled states in A4 for which two of the tangles 
{ t (ab)(cd),T(ac)(bd),T(ad)(bc)} equal to the maximal 
value of 3/2 while the remaining tangle equals to 1. 

Theorem 11. Let \ip) £ A, and denote by {Pj}, {Qj}, 
and {Rj} the eigenvalues of the reduced density matri- 
ces of \ip) obtained after tracing out qubits C and D, B 
and D, and B and C, respectively. If the Shannon en- 
tropies H({Pi}) = H ({Qi}) = 2, then H ({R t }) < I, 
with equality if and only if up to a local unitary \tp) — 
\C 2 )£M. 

Note that from the theorem above it follows that, up 
to a local unitary, \ip) = \C\) ( \ip) = IC3)) if one replace 
the condition H ({Pi}) = H ({QJ) = 2 with H ({Qi}) = 
H({Ik}) = 2(H({P i }) = H({R i })=2). 

Proof. Let \ip) £ A, and define two orthogonal 4x4 or- 
thogonal matrices: 




Theorem 12. Let ip £ T-L^ and a > 2. Then, 

Ei a) (4>) <5/3, 

with equality if and only if up to local unitaries, ip is one 
of the cluster states given in Eas.li3lj T5\) . 



Proof. We first prove it for the case a — 2. In this case, 
the Renyi entropy of degree 2 (which also called colli- 
sion entropy) can be expressed in terms of the linear en- 
tropy. This implies that the Renyi entropy of entangle- 
ment E^~ 2 ^ can be expressed in terms of the tangle: 



4 Q=2) (i^ s )) = -iog(i-^G^ B )) 



To simplify notations, we denote by t\, ti and t-$ the 
values of t { ab K cd)(^), T( A c)(bd)W, and T {AD){BC) (ip), 
respectively. With these notations we have 



E. 



(a=2) 



log 



Note that since T2 < 4/3 we have t\ + ti + tj, < 4. Now, 
the function — log(l — t/2) increase with t. Therefore, 
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E 2 2 obtains its maximum value when t\ + t% + £3 = 4. 
Now, we define 

on the domain 

D = {(ti,t a )|l <*i < 3/2 , 5/2 -ti < < 2 < 3/2}. 

A simple analysis of the function f{t\,t 2 ) implies that 
/(fiifo) obtains its minimum value of 1/32 only at the 
points (1,3/2), (3/2,1), and (3/2,3/2). Therefore, for 
states with these values for t\ t 2 and t 3 — 4 — t\ — t 2 , 
E 2 a ~ 2 ^ obtains its maximum value. In theorem [TT] we 
have seen that the 3 cluster states are the only ones with 
these values of t 1: t 2 , and t 3 . This completes the proof 
of the theorem for a = 2. The case for a > 2 follows im- 
mediately from the fact that the Renyi entropy is a non- 
increasing function of a and therefore E 2 a ~ 2 ^ > E 2 a>2 \ 
To complete the proof, we observe that for the cluster 
states E [ 2 a) = 5/3 for all a. □ 

V. MORE MAXIMALLY ENTANGLED 
4-QUBITS STATES 

In this section we characterize all the states in ip £ 
A that maximize the average a-entropy of entangle- 
ment E 2 a ^ (as defined in Eq. (IT9)) ). for given values of 
T(AB)(CD), T(ac)(bd), and T(ad)(bC)- We also character- 
ize the states in M. for which one of the three tangles 
t (ab)(cd), T(ac)(bd), and T(ab){bc), is equal to 3/2. 

We start with a definition of a class C of maximally 
entangled states: 

C = 

^e iS u + \f^J2 «j| l^P^l, cos 2 6 < I 

Proposition 13. Up to local unitaries, the class C con- 
sists of all the states in A such that all three distributions 
{Pj}, {Qj}> and {Rj}, as defined in Eq. \21$ , have the 
form {po>Pi)Pi)Pi} with po > pi; i.e. they are of the 
form given in Ea A2U\) . 

Remark. Note that the state \M) belongs to C. It corre- 
sponds to 9 — ir/2 and p — 1/2. 

Proof. Let ip — zqUq + Zitii + z 2 u 2 + Z3U3 be a state 
in A with the properties mentioned in the proposition. 
Therefore, since Pj = \zj\ 2 , up to local unitaries 

with p > 1/2. From the definitions of {Qj} and {Rj} in 
Eq. (|2~4"1) . and from the requirement that three of the QjS 



are equal and smaller than 1/6 and also three of the RjS 
are equal and smaller than 1/6, we get that 9\ — 9 2 = 
#3 = 6, and cos 2 9 < (l—p)/3p. Therefore, up to a global 
phase ip £C. □ 

Corollary 14. Let tp £ H4 and 4> £ C. If 

T (AB)(CD)(*P) < T(AB)(CD)W 
T (AC)(BD)W < T(AC){BD){4>) 
T (AD)(BC)W < T (AD)(BC)( < I ) ) , 

then 

E { 2 \tP) < E { 2 a) {0) for0<a<2 

E { 2 \^) > E { 2 \$) for a > 2 (32) 

with equalities if and only if <f> — tp up to local unitaries. 

The corollary follows directly from the proposition 
above and from theorem [5] 

Note that the corollary above also implies that for < 
a < 2, 

E< a \l>)<4 a \\M)), 

for all -0 € H 4 with T(ab)(cd)(^), t ( a C)(bd) (^), and 
t (AD)(bc) (VO all being no greater than 4/3. 

We end this section by classifying all the states in M. 
for which one of the tangles T(ab)(cd), t {ac)(bd), and 
T(ad)(bc) is equal to 3/2. Without loss of generality we 
will assume that tiab)(gd) — 3/2. 

Proposition 15. Let ip £ M. If T(ab)(cd) (VO = 3/2 
then, up to local unitaries, 

1 i e ie ie l ° 
H) = 2 U ° + 2 Ul + ~ U2 + ~ m ■ ( 33 ) 

Remark. The two cluster states with T(ab)(CD) =3/2 
have this form with 9 — Qot9 — tt/2. Among all the 
states of this form, the choice 9 = 7r/4 gives the high- 
est value for the average entropy of entanglement, but 
it does not reach the average entropy of entanglement of 
the state \M) (see Fig. 2). 

Proof. From Tiab)(gd) =3/2 we have 

3=0 

Now from the condition in Eq. (fTB"]) , we get Ylj=o e%28i = 
0. This condition with the freedom of global phase leads 
to the form in Eq.®. □ 
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' 0.5 1 1.59 2 2.5 3 



a 

FIG. 2: A graph of the average Renyi a-entropy of entangle- 
ment as a function of a. The blue line corresponds to the 
state \M), the green line to the state \L), the dashed red line 
to the cluster states, and the dashed purple line to the state 
given in Eq. (|33l) with 9 — 7r/4. For a > 2 we proved that the 
cluster states maximize However, we conjecture that 

the cluster states maximize E^" for a > Qo = 1.59... while 
the state \M) maximize it for a < Qo- Among all the states 
in M the state \L) and \M) minimize for < a < 2 and 
a > 2, respectively. 



VI. CONCLUSIONS 

Four qubits entanglement is far more complicated to 
analyse than its three qubits counterpart. This intricacy 
manifests itself with the uncountable number of inequiva- 
lent SLOCC classes. Such complexity also occur in 5 and 
6 qubits systems, although for these systems there exist 
maximally entangled states (such as the 5-qubits code 
state) with the property that any bipartite cut yields a 
maximally entangled (bipartite) state. Since such states 
do not exists in 4-qubits nor in n-qubits with n > 8, the 
study of 4-qubits entanglement gives an insight to the 
structure of n-qubits maximally entangled states with 
large n. Indeed, some of the results presented here, such 
as theorem O can be extended to n-qubits [33| . 

In this paper we found an operational interpretation 
for the 4-tangle as a kind of 4-party residual entangle- 
ment that can not be shared between 2-qubits/2-qubits 
bipartite cuts. This operational interpretation enabled 
us to find a family of maximally entangled states that 
is characterized by four real parameters. All the states 
in the family maximize the average bipartite tangle, but 
only two states in the family (i.e. the states \M) and \L) 
in Eqs. (|2ll[ )), maximize all the average Tsillas a-entropy 
of entanglement. In this sense, up to local unitary trans- 
formations, there are only two maximally entangled four 
qubits states. 

Both the states \M) and \L) are symmetric; that is, 
up to local unitaries, they are both invariant under per- 
mutations of the 4 qubits. The eigenvalues of their re- 
duced density matrices that obtained after tracing out 



two qubits have the form given in Theorem[8l Therefore, 
since they both maximize the average bipartite tangle, we 
believe they also optimize many other averages of bipar- 
tite entanglement monotones that were not introduced 
here. Moreover, the techniques introduced here suggest 
that states with the properties of_|M) and \L) may exists 
in higher dimensional systems [33j. 

We also found that the three cluster states in 
Eqs. (|3l4l5l) are the only (up to local unitaries) 4-qubits 
states that maximize the average tangle, and have the 
property that out of the three reduced density matrices, 
that obtained by tracing out two qubits, two are propor- 
tional to the identity. In addition, we showed that the 
cluster states optimize the average Renyi cv-entropy of 
entanglement with a > 2. The reason that it is the clus- 
ter states and not \M) or \L) that optimize this average 
Renyi entropy, is that the Renyi entropy with a > 2 is 
not concave and so the Renyie entropy of entanglement 
is only a deterministic entanglement monotone and not 
an ensemble monotone. 
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Appendix A: Kempf-Ness Theorem 

The purpose of this appendix is to state the version 
of the Kempf-Ness theorem that is used in this pa- 
per. Let H n be n-qubit space and let G be the sub- 
group SL(2, C) <8> SL(2, C) ® • • ■ ® SL(2, C) (n-copies) in 
GL(H n ). Let K = SU{2) ® SU(2) <g> • • • <g> SU{2). Let a 
be the Lie algebra of G contained in End(H n ). We set 
Crit{H n ) = {0 G H n \ (4>\X\4>) = 0,X G Q}. The Kempf- 
Ness theorem in this context says (the only hard part of 
the theorem is the "if part of 3. which we don't use in 
this paper): 

Theorem 16. Let <j> G H n then 

1. (f> G Crit{H n ),g G G then ||<7</>|| > ||<^|| with equality 
if and only if gcj> G K <j>. 

2. If <p G Hn then (f> G Crit(H n ) if and only if \\g4>\\ > 
\\<f>\\ for all g EG. 

3. If 4> G "Hn then G(f> is closed in T-L n if and only if 
G(l>nCrit{H n ) 5^0. 

We now assume n = 4. Let A be as in section II. Then 

Proposition 17. Crit(H 4 ) = KA. 

Appendix B: Proof of theorem [8] 

We first prove the theorem for the von-Neumann en- 
tropy (i.e. the case a = 1) and then we will consider 
the case a ^ 1 separately. We want to optimize the 
function /(A , Ai, A 2 , A 3 ) = - J2l=o ^ log A fe under the 
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constraints J2k=o Afc = 1 an d J2k=o A fc = 1 — r /2> while 
< Afc < 1. The Lagrangian is therefore given by 



3 / 3 

£ = -^A fc logA fc +,i( VA fc I 

k=0 



\fe=0 / 



where /i and are the Lagrange multipliers. Therefore, 
the critical points in the interior of the domain (i.e. < 
Xk < 1) must satisfies the equation: 

dC 

— = -logA fc -loge + /i + 2//A fe = (Bl) 

Now, we first show that if all for Afc satisfy the equation 
above, then the set {Afc} contains at most two distinct 
numbers. To see that, suppose that there are three dis- 
tinct numbers. Then, without loss of generality, lets as- 
sume that Ao > Ai > A2 > 0. Thus, from the 3 equations 
above (for k = 0, 1, 2) it follows that 

(A - Ai) log A 2 + (Ai - A 2 ) log A = (A - A 2 ) log X 1 

Denote by a = (Ai — A 2 )/(Ao — Ai). Hence, a > and 

log A 2 + a log Ao = (1 + a) log Ai , 

which is equivalent to 

\ \a \l+a 

Denote by x = A 2 /Ai and y = A0/A1. Therefore, x < 1, 
y > 1, a = (1 - x)/(y - 1), and 

xy a = 1. 

From the last equation and the generalized arithmetic- 
geometric mean inequality we get 

1 = (x^) 1 /^) < -L_ (x + ay) = 1 , 
1 + a 

where the last equality is obtained by substituting a — 
(1 — x)/(y — 1). Not that the geometric-arithmetic mean 
inequality is saturated if and only if x = y and therefore 
we get a contradiction to the assumption that there are 
3 distinct numbers in the set {Ai}. Therefore, for the 
interior points we have 3 options: (a) Ao > Ai = A 2 = 
A 3 , (b) A = Ai > A 2 = A 3 and (c) A = Ai = A 2 > 
A3. Option (b) is the only one that does not appear 
in theorem [HI In this case r must be greater than 1, 

and A 2 = A 3 = K 



and A = Ai = 1+V ^ (r) 
is a straightforward calculation to show that the von- 
Neumann entropy of this distribution never equals the 
von- Neumann entropy of the distributions that appear in 
the theorem. (Note that this is all we have to show since 
it has already been proved in [3l| that the distributions 
in the theorem are the optimal ones). 

As for the critical points on the boundary, set A3 = 
and then the same argument as above implies that the 
set {Aq,Ai,A 2 } contains at most two distinct numbers. 



We therefore have two options: (a) A = Ai > A 2 and 
(b) A > Ai = A 2 . Again, the distribution (b) does not 
appear in the theorem, but it is straightforward to show 
that its von-Neumann entropy never equals to the von- 
Neumann entropies of the distributions in the theorem. 
The last point on the bounday that we need to check is 
when A3 = A 2 = 0, but this distribution appears in the 
theorem. This completes the proof for the case a = 1. 

We now prove the theorem for the case a 7^ 1 (as 
well as a ^ 2). In this case, we optimize the function 
/(Ao, Ai, A 2 , A3) = Y^k=o^k un der the same constraints 

above; that is, J2k=o^ k ~ 1 anc ^ J2k=o ^fe = 1 ~ r /^- 
The Lagrangian in this case is given by 



k=0 



\k=0 



\k=Q 



where \i and v are the Lagrange multipliers. The critical 
points in the interior of the domain must satisfies the 
equation: 



9C ! 

dTk = aXk 



loge + ^ + 2^A fe = 



(B2) 



Similarly to the argument above, we first show that the 
set {Xk} contains at most two distinct numbers. To see 
that, suppose that there are three distinct numbers Ao > 
Ai > A 2 > 0. Thus, from the 3 equations above (for 
k = 0,1,2) it follows that 



A o 



A" 1 Ao — Ai 



\ot— 1 \ot- 
A ~ A 2 



Ao — A 2 



Now, denote by x = A1/A0 and y = A 2 /Ao. From our 
assumptions < y < x < 1. With these notations the 
equation above can be written as: 



1 ~x c 



i-y° 



l-x 



1-1/ 



However, for non- negative a 2 the function f(x) — 
(1 — x Q_1 )/(l — x) is one-to-one and therefore we get 
a contradiction. This complete the proof that the set 
{Afe} contains at most two distinct numbers. The rest of 
the proof follows the same lines as the proof for the case 
a = 1. 



Appendix C: Global max of /, 

J min 



(a) 
max 



and global min of 



1. Calculation of U (a) 

In this section we prove that for a > 2, U^ a ' = 
/max(4/3,4/3,4/3), and (4/3,4/3,4/3) is the only point 

of global maximum. Since the function /max(ii, fe, is) is 
invariant under permutations of t\, t 2 , and £3, we can as- 
sume without loss of generality that the global maximum 
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of /max is obtained at a point with h > t 2 > t 3 . We will 
therefore look for the maximum of /maL(ii > *2> 4 — ti — ia) 
in the domain 



D 



(*1,*2) 



4 3 „ 

- < ti < - , 2 
3 ~ 1 - 2 ' 



1 



-ti <u<u 



In this domain, £2 can be either bigger or smaller than 
4/3 and therefore the derivative of /max with respect to t% 
is not continuous at points with t 2 — 4/3. We therefore 
split the domain D into two regions 



4 3 4 

-<*].<-, - <t 2 <t 1 
1 , _ 



£»2 = ((ti,t 2 : 



4 3 „ 1 4 

3^*^2' 2 -2*^*^3 



so that on D\ (or D 2 ) all the derivatives of /max are 
continuous. We start by maximizing /max on the domain 
Di. 



a. Maximizing /,£L on t/ie domain D\ 



Denote by Xi 



1 — %ti for z 



1,2 and y 



1 — §£3. In these variables, the function gmaxO^i , £ 2 ) 



frnJxih, t 2 , 4 - £ x - i 2 ) is given by 



7 («) 

Ml.! ax 



(xi,x 2 ) 



a- 1 



E 

i=l,2 



(l + a; 4 ) Q , 1(1-3^)" 



4" 



- 3 ^( 1 +w) a - 3 ^( 1 - 2 ») tt 

where in term of the variables x\ and x 2 , y — 
I y^l — + In terms of these new variables, the 

domain Z?i is given by < x\ < x 2 < 1/3. We start by 
looking at the critical points in the interior of D± . 

The critical points of gmlx(xi,x 2 ) satisfies the condi- 
tions 



3 fJlXltlX ' i 

-(xi,x 2 ) = — (u a (y) - v a (xx)) = 



dx\ 



Q ™^ {xi,x 2 ) = — (u a (y) - v a (x 2 )) = , 



where 



Hence, the point (xi, x 2 ) is critical if and only if v a (xi) = 
v a (x 2 ) = u a (y), where y = j^/l — + x\) (note that 

< y < 1/2). In the next two lemmas we prove two 
useful properties of the functions v a (x) and u a (y). 



Lemma 18. If 2 < a < 4, and < y < 1/2, t/ie« 

<(y) < 0. // 2 < a < 5, and < x < 1/3, then 
v' a (x)<0. 

Proof. A simple calculation gives 

u'M = 

(a-^.-i ^ t (1 + - 1} + (1 - 2y) " (1 + ^ ' 

where /3 = a — 2. From the assumption of the lemma < 
/3 < 2. Therefore, one can easily check that u' a (0) < 
and u^(l/2) < 0. All that is left to show is that in the 
domain < y < 1/2 the function 

w{y) = (1 + yf{py - 1) + (1 - 2yf(l + 2(3y) 

is negative. To find its maximum value, we calculate 
its critical points. The requirement w'(y) = gives (1 + 
yf- 1 = 4(l-2y)' 3 - 1 . Clearly, there are no critical points 
for P < 1 in the domain (0, 1/2). For f3 > 1 we express 
the value of w(y c ) at the critical point by substituting for 
(1 + 2/ c ) /3 ~ 1 the value 4(1 - 2y c ) /3 ~ 1 . This gives, 

w(y c ) = -3(1 - 2y c f- 1 [1 - 2{[3 - l)y e ] < , 

for y c < 1/2 and (3 < 2. Hence, u' a (y) < for 2 < a < 
4. Following the same arguments, one can show that 
v' a {x) < for 2 < a < 5. □ 

Lemma 19. If a > 4 i/ien f/ie global maximum of v a (x) 
(in the domain < x < 1/3) is strictly smaller than 
the global minimum of the function u a {y) (in the domain 

< z/ < 1/2;. 

Proof. The global extremum points of u a and v a are ob- 
tained on the boundary or on critical points. Therefore, 
we first check the bounday (that is, end points). We have 



«a(0) 



o 



4«-2 

a 



a 



1 



a - 1 3 Q - 2 
a 1 



Clearly, for a > 4, we have max{u a (0), v a (l/3)} < 
min{?/ Q (0), it Q (l/2)}. We now estimate the values of u a 
and v a at their critical points x c and y c . From u' a (y c ) = 
and v' a (x c ) = we have 



«a(lfc) = 3^1 [(1 + ^c)"- 2 + 2(1 - 2y c y 



««W = t^t [(i + ^) Q " 2 + 3(i - 3x c ) 



(Cl) 



Since we do not have explicit expressions for x c and y c , 
we find an upper bound for v a (x c ) and a lower bound for 
u a (y c ). Since the functions in Eq. (|Clj) are convex for 
a > 4, we get 

v a {x c ) < max{u a (i c = 0), v a (x c = 1/3)} 
la a 1 



3 a-2 ' 4 a-2, 
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The minimum value of the function u a (y c ) given in 
Eq. (|C1I) is obtained at the point 



2/c = 



41/(0-3) _ i 
1 + 2 • 4 1 /("-3) 



Note that this is not necessarily the true value of y c , 
but rather the value at which the function in Eq. (|Cip is 
minimized. It is a straightforward calculation to show 
that at this value of y c 

u a {y c ) > max j-^^ , v a (0) , v a (1/3)} • 

This completes the proof that u a (y) > v a (x) for a > 4 
and for all x G [0, 1/3] and y G [0, 1/2]. □ 

From lemma [19] it follows that for a > 4 the function 
gmax(%i, x 2 ) does not have critical points in the interior 
of D\. For 2 < a < 4, gmax{xi,%2) can have critical 
points. However, from the second derivatives test and 
from lemma [IS] it follows that the Hessian is positive 
definite. Therefore, these critical points are local min 
and can not be a global max. The global maximum of 

(/max(xi , 

x 2 ) is therefore obtained at the boundary of D\. 
The boundary of D\ is a triangle with 3 sides given by 
x\ = 0, Xi = x 2 , and X2 — 1/3. If Xi = then 



dgmax(0,X 2 ) 
» 



X 2 



(u a (y) - v a {x 2 )) 



Therefore, <7ma!x(0, X2) is convex for 2 < a < 4 (see 
lemma 1181) . and has no critical points for a > 4 (see 
lemma [H]). Therefore, its global max is obtained at one 
of the end points (0,0) or (0, 1/3). 
On the side x± = x 2 = x we have 



( x > x ) = 77 (u a (y) - v a (x)) 



Hence, the same arguments implies that the global max 
of <?ma X (x, x) is obtained at one of the end points x = 
or x = 1/3. Similarly, on the side x 2 = 1/3, the function 

ffmax(xi, 1/3) obtains its global maximum at one of the 
end points x± = or x\ = 1/3. Among the three vertices 
(0,0), (0,1/3), (1/3,1,3), we have 

3^1(1/3, 1/3) > max {^1(0,0) , ^(0, 1/3)} 

for all a > 2. Hence, on D\, g^lx obtains its global max 
at the point (xi,Xa) = (1/3, 1/3) which is equivalent to 
h = t 2 = t 3 = 4/3. 



b. Maximizing /ril^c on the domain D2 



Denote by x = \/l — |ti and yt = \/l — jti for 



2,3 . In these variables, the function hmax(lf2, Va) 



/i a ax(4 - f 2 - *3, t 2 , t 3 ) is given by 



^mix(2/2,2/ 3 ) = 

1 



1 



3 «+i 

i=2,3 



(l + xr + -(l-3x) a 



where in term of the variables y 2 and j/3, x = 
— 8(2/2 + 2/1)- l n terms of these new variables, the 



domain D 2 is therefore given by 



{(2/2, 2/3) y 2 + vl < g > 2/3 > 2/2 > 0} . 



£ 2 = 



We first look at the critical points in the interior of D 2 . 

The critical points of /imax(j/2, 2/3 ) satisfies the condi- 
tions 



<9?/2 
#2/3 



(2/2,2/3) 



(2/2,2/3) 



2?/2 



2ys 
9 



(i>„(x) - u a {y 2 )) = 



(l> a (x) - W Q (j/ 3 )) = 



Hence, it follows from lemma IT51 and lemma [1^] that the 
global max of hmlx(y 2 , 2/3) is obtained on the boundary 
of D 2 . 

The boundary of D 2 consists of the line y 2 = 0, the line 
2/i = 2/2, and the curve (2/2,2/3) = (sm9/V&, cosd/yS) 
with < 9 < tt/4. The global maximum of hmax on 
the lines 2/2 = and j/i = y 2 is obtained on one of the 
endpoints (0,1/2), (0,0), and (1/4,1/4). The argument 
follows from lammas ll8l and [T9l in the same way as it was 
used in the analysis of the boundary of D\. Therefore, 
we focus now on the curve (2/2, 2/3) = (sin 0/\/8, cos9/y/E) 
with < 9 < tt/4. 

Let 



HO) = h<£l 



sin 9 cos 9 



Note that for these values of y 2 and 2/3, x = 0. Hence, 
sin(20) 



72 



cos f 

7f 



sm f 

71" 



From lemma [15] the function u a is one-to-one for 2 < 
a < 4. Therefore, the only critical point in this case is 
(2/2,2/3) = (1/4, 1/4). For a > 4 it is a simple calculation 
to verify that h{9) < htL(0,0). Therefore, since 

^2x(0,0) > max {^(1/4, 1/4) , ^ x (0,l/2)} 

for all a > 2, we conclude that on D 2 , hmax obtains its 
global max at the point (2/1,2/2) = (0,0) which is equiva- 
lent to h=t 2 = t 3 = 4/3. 
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2. Calculation of L (a) 

In this section we prove that for < a < 2, = 
.^(4/3,4/3,4/3), and (4/3,4/3,4/3) is the only point 
of global minimum. From theorem [5] it follows that 
/min(*i'* 2 '*3) for < a < 2 is given by the exact same 
expression as /nSL(*i, *2j £3) f° r a > 2. Therefore, our 
proof that L( Q ) = 7^(4/3,4/3,4/3) follows the exact 
same steps used in the calculation of U^ a > for a > 2. 
The only difference is that the lemmas [18] and [19] do not 
hold for < a < 2, and instead we have the following 
lemma. 

Lemma 20. If < a < 2, < x < 1/3, and < y < 
1/2, then u' a {y) > and v' a (x) > 0. 

Proof. A simple calculation gives 



1 



(a- 



1)3 Q - 1 y 2 



(l + 20y) (1 + Pv) 
(l-2y)P {l + yf 



where ft = 2 — a. From the assumption of the lemma 
< j3 < 2. Therefore, one can easily check that 
linij,_j.o u' a (y) > and lim a _>.x/2 u' a (y) = +00. All that 
is left to show is that in the domain < y < 1/2 the 
function 



w{y) 



I 



1 



(1 - P) y 2 



(1 + 2j3y) (1 + py) 
(l-2y)P (1+yf 



is positive. To find its minimum value, we would like 
to calculate its critical points. However, the requirement 
w'(y) = gives 4(1 + yY +1 = (1 - 2y) p+1 . Hence, there 
are no critical points in the domain (0,1/2). That is,, 
u' a (y) > for < a < 2. Following the same arguments, 
one can show that v' a (x) > for < a < 2. □ 

With this lemma replacing lemmas [T51and[T9l the proof 
that = /^(4/3,4/3,4/3) follows exactly the same 
steps that appear in the calculation of U^ a \ 
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Note that /o is the G-concurrence [231 ] between qubits 

(1,2) and (3,4) 

A state ip is stable if the set of states Gip is closed in H4. 
We are willing to conjecture that it is true for all a > 0. 



